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We present general solutions to the equations of motion for a superconducting relativistic chiral 
string that satisfy the unit magnitude constraint in terms of products of rotations. From this result 
we show how to construct a general family of odd harmonic superconducting chiral loops. We further 
generalise the product of rotations to an arbitrary number of dimensions. 

PACS numbers: 11.27.+d,98.80.Cq 



I. PRELIMINARIES 

Particle physics models where symmetry breaking is 
involved predict, in many cases, the existence of topolog- 
ical defects, which are formed when the topology of the 
vacuum manifold of the low energy theory is non-trivial 
jij. Cosmic strings, in particular, are line-like objects 
that are formed when the the vacuum manifold contains 
unshrinkable loops. For a review see 

In H it was shown that cosmic strings can be super- 
conducting. In the case when the charge carriers on the 
string are not coupled to a gauge field the action for the 
string and the current can be taken to be 

s = J d 2 e^(-A* + li ab M b ) (i) 

where /i is the mass per unit length of the string, 7 ab is 
the induced metric on the string worldsheet and <f> is the 
field of the charge carriers living on the string. These 
strings were shown in 0] and || to have solutions in the 
case when r ) ah <t> ,a<t> ,b — of the form 

x=i[a(u) + b(«)] (2) 
for the string position and 

<t>=\f(v) (3) 
for the field living on the string with the constraints 

a' 2 = 1 (4) 

and 

b' 2 + /' 2 = 1 (5) 

where u = a — t and v = a + r and a and r are space-like 
and time-like parameters respectively that parametrise 
the string world-sheet. These strings are called chiral 
because the current only moves in one direction on the 
string. 



Comparing this to the usual Nambu-Goto case one can 
see that f(v) acts like a fourth component of the three- 
vector b(w), making chiral superconducting strings be- 
have like Nambu-Goto ones with chiral excitations in an 
extra fifth dimension. Indeed, this property was used in 
in an investigation of the properties of superconduct- 
ing cosmic string cusps. 

The right- and left-moving excitations, a' and b', on a 
regular Nambu-Goto string in Minkowski space-time are 
arbitrary functions that satisfy the unit magnitude con- 
straint, | a' | = |b'| = 1. Expressions for these functions 
are often given as Fourier sums and the unit magnitude 
constraint generally gives a non-linear set of equations 
involving the vector coefficients of the Fourier expansion. 
As a result, parametrising strings beyond the first few 
harmonics proves to be a difficult task. Fortunately, in 
that case, there exists a method to generate strings in- 
volving products of rotation matrices Jjj that act on a 
starting unit vector so that the unit magnitude constraint 
is satisfied trivially. 

In a recent study of the properties of chiral cosmic 
strings Q it was assumed that the current is constant. 
The work in ]9| assumed that the current takes a very 
simple non-constant form. As was pointed out in the 
latter work one could expect to have loops with varying 
currents if the loops are formed by intersections involv- 
ing different strings or if different segments of the loop 
or string were at some point in causally disconnected re- 
gions. For long strings this is always the case and we 
therefore expect varying chiral currents to be generic. 

The purpose of this work is to generalise the work in 
[Q for generating three dimensional unit vectors to four 
dimensional ones that include the current as a fourth 
component of b(cr + r). We start by casting the method 
somewhat differently and generalise it to four dimensions. 
From this result we construct a family of chiral supercon- 
ducting odd-harmonic loops. We further generalise the 
product of rotations to an arbitrary number of dimen- 
sions. 

In Section 2 we show how to construct an arbitrary 
N harmonic unit vector in four and three dimensions. In 
Section 3 we use these results to construct arbitrary chiral 
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current carrying superconducting odd harmonic loops. 
In Section 4 we generalise the arguments in Section 2 to 
an arbitrary number of dimensions and we conclude in 
Section 5. 



II. SOLUTION TO THE UNIT MAGNITUDE 
CONSTRAINT IN TERMS OF PRODUCTS OF 
ROTATIONS 



An ■ Aat_i — Bjy ■ Bjv-i — 0, 
An ■ Bat_i + • Aat_i = 



(13) 



respectively. It follows from (12) that the highest har- 



monic is a circle that lives on some arbitrary plane. 
Clearly, we can introduce coordinates such that An and 
B j\r lie on the w and x axes 



aw, B 



N 



(14) 



We can think of the Euclidean 4-vector 



b' 



K 
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making the highest harmonic a circle of radius a on the w- 
x plane. This puts the vector b^y into the so-called stan- 
dard form |?]]. Let R wx {9) be a matrix that rotates the 
^ w-x plane by an angle 9. Acting on b^, with, R wx (— v) 
in these coordinates lowers the highest harmonic term of 



as having unit magnitude according to (Jfy with b w = f. 
Consider the vector b^, that can be constructed from a 
finite sum of Fourier components, 



N 



o N (v) = Z + y^{A n cosnw + B„ sinra;}. 



n=l 



(7) 



The Fourier coefficients satisfy the set of 4TV + 1 non- 
linear relations derived in |7| 



JV 



E ( 



n—m — N 

with m = 0, 1, . 



(X n • OLr 



,2N, 



AS m0 



(8) 



N 

E 

n=m—N 



Pn ■ a m - n ) = 



(9) 



with m= 1, ...,27V, 



a- n = A^, p n = -/3_„ = B„, n / 0, (10) 



and 



a = 2Z, /3 = 0. 



(11) 



These equations can be obtained from the constraint 
equation (||). The total number of degrees of freedom 
in the coefficients in (Q) is 8/V + 4 so the remaining num- 
ber of degrees of freedom after satisfying the constraint 
is AN + 3. Below we show how to construct (0) from a 
product of rotation matrices by generalising a modified 
version of the three dimensional method presented in 
to four dimensions. 

The constraint equations (||) and (||) for m = 2N and 
m — 2N — 1 are 



A% — B>v, An ■ Bjy — 



(12) 



and 



Rwx(-V) 



( acosNv 
a sin Nv 


V 



/ a cos(7V — l)v 
a sin(7V — l)v 


V 



(15) 



This is not sufficient to verify that the overall harmonic 
content has been lowered because the N — 1 terms of (Q) 
could still give us an N harmonic through trigonometric 
identities. We now show, however, that this is not the 
case. 

The constraint equations for m — 27V — 1 ( |l3|) give us 
using ( |l4| ) the conditions on the coefficients, 

(Ajv-i)ui = (Bjv-i)k, (An~i) x — — (Bjv-i) w . (16) 



It is not too difficult to see that using the conditions dl£ ) 
when acting with R wx (—v) on the TV — 1 terms of (|7) 
does not lead to the creation of TV harmonic terms, 



R wx (-v) [A N -! cos(TV - l)v + B N -i sin(7V - l)v] (17) 

(A N -i) w cos(TV - 2)v - {A N -i) x sin(TV - 2)v 
{A N -i) x cos(TV - 2)v + {A N -i) w sin(TV - 2)v 
(A N -i) y cos(TV - l)v + (Bjv-i) v sin(TV - l)v 
{A N -i) z cos(TV - l)v + (B A r_ 1 ) 2 sin(7V - l)v 

Clearly, the TV — 2 terms will not yield an TV harmonic 
term when acted on by R wx {—v). 

We now have a string in the form of (0) with TV — * TV — 
1. Its highest harmonic is therefore also a circle of some 
other radius living on some other arbitrary plane. Armed 
with this knowledge we can see that the unit magnitude 
constraints are such that in general we can write 



V n = Rp n (v)V } 



N-l 



(18) 



where Rp N (v) is a rotation by an angle v on the plane 
Pn where the highest harmonic lives. By induction it 
must be that 



b' N =l[Rp t (v)b' 



(19) 
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where the Rp i (v) are rotations by an angle v on arbitrary 
planes Pi and h' is an arbitrary constant unit vector in 
4 dimensions. 

To specify a plane in 4 dimensions one needs to specify 
a direction on the plane (3 angles), a linearly independent 
direction (2 angles) , but now the plane is overspecified by 
internal rotations (—1 angle), giving a total of 4 degrees 
of freedom for the matrices Rp i (v). For an N harmonic 
vector, one therefore has 4iV parameters in the rotators 
and 3 parameters in the constant unit vector h' giving a 
total of AN + 3 parameters which checks perfectly with 
8iV + 4 vector coefficients in (Q) minus the AN + 1 con- 
straints (§) and @. 

The form of the Rp i {v) matrices is quite simple. Gen- 
erally, if we want to transform, say, a rotation by an angle 
v on the w-x plane to one on an arbitrary plane we will 
need to perform a transformation of the type 



R Pi (v) = E t R wx (v)E? 



(20) 



For the purpose of finding the form of the rotators Ei it 
is easiest to envision the inverse process to the one we 
are seeking, namely the rotation of an arbitrary oriented 
plane to lie on the w-x plane. Let's consider first the 
projection of the four dimensional arbitrary plane onto 
the (x,y,z) subspace. We can perform a rotation by 
some angle a about the z axis (R xy (a)) until the vec- 
tor perpendicular to the projected plane lies on the y-z 
plane and perform a further rotation by an angle f3 about 
the x axis (R yz (f3)) until that vector lies on the z axis. 
At this stage the projected plane lies wholly on the x-y 
plane. The ranges of both a and j3 from to tt are suffi- 
cient to perform these transformations. After performing 
these two rotations, our original four dimensional plane 
lies entirely in the (w, x, y) subspace and we can repeat 
an analogous process to the one above to rotate it into 
the w-x plane. A rotation by an angle 7 about the y 
axis (R W x("f)) puts the vector perpendicular to the plane 
on the x-y plane and a rotation by an angle S about w 
(R xy (S)) makes that vector parallel with the y axis. For 
the first of these rotations a range of 7 from to tt is 
sufficient, for the second rotation, however, matters are 
slightly different. If the plane we were trying to rotate 
was featureless it would be enough for the range of S to 
be from to tt. In fact this is not the case. The plane 
contains a circle in v which can be oriented clockwise or 
anti-clockwise on the w-x plane and therefore the final 
rotation on the x-y plane in general requires an angle 
that ranges from to 2tt. 

Keeping in mind these considerations we can quite gen- 
erally write 



Ei — Rxy{ai)Ryz{Pi)Rwx(li)Rxy{5i) 



(21) 



where on, 74 range from to 7r and Si ranges from 
to 2tt. Then 



bjv(«) = II E l R wx {v)Ejh' . 



(22) 



In order to construct the entire chiral string, we also 
need to find the form of a(u) in (|J). The constraints 
([I|) can be satisfied using a product of rotations that can 
be found from analogous arguments to the ones in the 
preceding section. For M harmonics this yields 



ak(«) = II D t Rxy(u)Df a ' 



i=M 



with the rotator 



Di — R x y 



)R 



yz \ 



i) 



(23) 



(24) 



where the angles 9i range from to tt, the angles fa from 
to 2ir and R are the three dimensional rotation matrices. 



III. AN APPLICATION: A FAMILY OF ODD 
HARMONIC SUPERCONDUCTING CHIRAL 
LOOPS 

A. Overall Orientation Freedom 

Both expressions for the oppositely moving excita- 
tions on the string (22) and ( p3| ) include overall orien- 
tation freedom. In some applications, for instance self- 



intersection or gravitational radiation analyses, only the 
shape of a loop, and not its orientation, is important. In 
this case the inclusion of overall orientation freedom of 
the right-moving and left-moving excitations separately 
is unnecessary: All that matters is the relative orienta- 
tion between a' and the spatial part of b' . 

Overall orientation freedom of the loop is set by an 
Euler matrix Q that acts only on a' and the spatial com- 
ponents of b' and contains three angles. We can use 
this freedom to standardise the vectors in some way. We 
will choose the plane of the last rotation in a' and b' to 
contain some coordinate axis (this will enable us in the 
following sub-section to use the arbitrariness in the origin 
of u and v to eliminate more parameters) . 

The four-dimensional circle that constitutes the high- 
est harmonic of b' projected onto three-dimensional 
space typically looks like an ellipse. We can perform or- 
dinary three-dimensional rotations on it to put it in some 
convenient form. In particular, we can rotate the ellipse 
on the x-y plane until its major axis, say, lies on the y-z 
plane and further rotate on the y-z plane until it lies en- 
tirely on the z axis. We still have one more rotation left 
in Q which we choose to be a rotation on the x-y plane. 
This leaves the major axis of the ellipse on the z-axis but 
orients the circle of the highest harmonic in a' such that 
it contains the rr-axis. 

We now construct a four-dimensional planar rotator 
that contains the z-axis but is otherwise arbitrary. Start- 
ing with R wz (v), we see that we can perform an arbi- 
trary planar rotation on it in the (w, x, y) subspace that 
preserves the z-axis. This requires only two rotations. 



i—N 
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Therefore, to rotate about a plane that contains the z- 
axis but is otherwise arbitrary we can use 

R X y(7)R W x(S)R WZ (v)R wx (S)R xy {-j) (25) 

where 5 ranges from to it and 7 from to 2tt. This 
procedure eliminates two of the parameters in the last 
planar rotation of b'. 

In three dimensions to rotate by an angle u on a plane 
that contains the x-axis, but is otherwise arbitrary, we 
may use 



R yz {9)R xy (u)R yz (—9) 



(26) 



where 9 ranges from to 2tt. This procedure eliminates 
one of the parameters in the last planar rotation of a'. 
This leaves a' and b' in the form 

1 

b' N (v) = R» Pn (v) J] E t R wx (v)Efb' Q . (27) 

i=JV-l 

with 

R b PN (v) =R wx {-i)R xy {8)R wz {v)R xy {-8)R wx {- 1 ) (28) 
and 

1 

a M ( u ) — Rp M ( u ) J] D t R xy (u)Df a ' (29) 



i=M-l 



with 



R Pm {u) = R yz (9)R xy (u)R yz (-9). 



(30) 



B. The Origin of u and v 

The conditions on the coefficients of the highest har- 
monics (|l2| ) only specify the planar rotation up to a phase 
so that generally we can take R b N (v + j3) and R% f (u + a) 



in (|27D and (|29|). 

We consider the action of this extra planar rotation 
matrix on W N (v). It can be verified to be 



b' N (v) = R P Jv) [] R P ,(v)R b Pff (0)b' o . 



i=N-l 



where 



R P ;(v) = R b P ((3)R Pi (v)R b P (-P) 



(31) 



(32) 



with analogous expressions for &' M (u) . It is important to 
note that the effect of replacing R Pi (v) with R P >(v) is 
to make the same transformation on each of the planes 
of rotation, in other words, to rotate on some other set 
of planes. Since we can express any rotation on a plane 
using ( |2~l| ) however, the effect of the matrices Rp N ((3) and 
R P (a) on the rotators can be ignored. 



Since the planar rotations R b P]v and Rp M both contain 
one of the coordinate axes we can choose $ and a so that 
the z component of b' and the x component a vanish. 
This leaves them in the form 



cos V b 



Rp N iP)b' a = I ^IZol I =R xy (MRUe b )w (33) 







and 



R a p M (aK= cos9 a =R yz {9 a )y 



(34) 



where <pb and 9 a range from to 2it and 9^ ranges from 
to 7r. This means we can write a' and b' as 



*' M (u) = Rpju) JJ R Pl (u)R yz {9 a ) 



(35) 



-M-l 



and 



b' N (v) = Rp N {v) Yl R Pi (v)R xy (<t) b )R wx {9 b )w. (36) 



=JV-1 



C. The Center of Mass Constraint 

In order to construct string loops, apart from solving 
the unit magnitude constraint, we need to satisfy the 
center of mass constraint, namely that the loop should 
be closed and that we want to work in the rest frame of 
the loop. These constraints imply that the center of mass 
term must vanish, Z = in (0), for both a' and b'. In 
general this is an intractable problem because the center 
of mass term Z is a non-linear function of the angles 
in the rotation matrices. We can, however, solve this 
problem in an analogous way to that introduced in [ fl0| . 
The guiding principle behind such a construction is to set 
the starting center of mass terms to zero, and to apply 
further rotation matrices in such a way as to ensure that 
trigonometric identities never lead to the production of 
a zero harmonic. In practice this was done by using odd 
harmonics only and choosing some parameters so that 
the starting center of mass term is zero. Here we will 
proceed along similar lines. 

Generally, if the starting unit vector lies somewhere on 
the plane of the first rotation, i = 1 in (|3f]) and ([56]) , no 
center of mass term will be generated. Since our starting 
unit vectors are given by (33) and ([34j) we need to choose 
the first planes of rotation appropriately. 

Earlier we established that if we want a rotation by an 
angle u on a plane that contains the x-axis we can use 
fl26| ) . If we want an arbitrary rotation that contains the 
y-axis instead we may use 
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R x MR xy (u)R xz {-<f>). (37) 

If we now decide we want a rotation by u on a plane that 
contains a vector lying somewhere on the y-z plane, it is 
sufficient rotate this last matrix (|37]) on the y-z plane by 
whatever angle the vector makes with the y-axis. This 
is precisely the situation in (|34|). We therefore write the 
first rotation, i = 1 in ( p5| ) as 

R Pl (u) = R yz {e a )R xz {^ l )R xy {u)R xz {-^ l )R yz {-e a ) 

(38) 

where <f>\ ranges from to 2tt. 

To find the first rotator in (|3^) we proceed analogously. 
We want to find a rotation by an angle v on a plane 
that contains a vector lying somewhere in the (w, x, y) 
subspace, but is otherwise arbitrary, because this is the 
situation of ( |33| ) . If we start with a rotator that contains 
the w-axis but is otherwise arbitrary, taking say, 

R yz {li)R X y{5 l )R wx {v)R xy {~8 1 )R yz {- 11 ), (39) 

where 5\ ranges from to 7r and 71 ranges from to 27r, 
it is not hard to see that the first rotator in (^5|) must be 

Rp 1 {v) = R xy (8b)Rwx((f>b)Rv Z (ii)R X y{Si)Rw X {v) 

x R xy (-5 1 )R yz (-j 1 )R wx (-(f> b )R xy (-e b ). (40) 

Further rotations should be by 2u and 2v to avoid the 
production of center of mass terms through trigonometric 
identities. These considerations yield 

2 

b' 2N _ 1 (v)=Rp N (2v) J] E t R wx (2v)E? 

i=N-l 

x R X y(0 b )R wx (<l)b)R xz (ii)R X y(5i)Rw X {v)w (41) 

and 

2 

= R a p M ( 2u ) II DiR xv {2u)Dj 

i=M 

x Rxyie^RyMRxyi^y. (42) 



IV. GENERALISATION TO HIGHER 
DIMENSIONS 

Here we show the straightforward generalisation of our 
four dimensional argument in Section 2 to a unit vector 
living in arbitrary dimensions. The constraint equations 
( |) and @ as well as the argument leading to ( fl9|) and 
(20|) are independent of the number of dimensions the 
vector lives in. The only thing that changes with the 
number of dimensions is the parametrisation of E, the 
rotator that takes an arbitrary oriented plane to the plane 
of the first two coordinates. In the following, we label our 
spatial coordinates by the numbers 1 through d. 



If we consider the projection of the plane onto the sub- 
space given by the last three coordinates (d — 2, d — 1, d) 
one can see that a rotation by a d about the d axis 
(Rd-2,d-i(&d)) until the vector perpendicular to the pro- 
jected plane lies in the d-d— 1 plane followed by a rotation 
by an angle (3d about the d — 2 axis (Rd-i,d(Pdj) is suf- 
ficient to rotate the projected plane out of the d axis. 
To perform these transformations it is sufficient for the 
angles to range from to ir. We can repeat this proce- 
dure by moving up d — 2 times in the coordinates until 
the plane lies entirely in the 1-2 plane as desired, ensur- 
ing that the range of the angle in the very last rotation 
R2,3(@3) is from to 2n to account for the fact that we 
are dealing with an oriented plane. We can then write 
the rotator as 

E = Rd-2, d -i(a d )R d - 1 , d {P d )...R 1 ,2(a 3 )R 2 , 3 {P3)- (43) 

The number of parameters introduced by such a product 
is 2(d— 2) per harmonic plus (d— 1) parameters to specify 
the initial unit vector giving a total of 2N(d — 2) + d — 
1 which checks with d(2N +1) degrees of freedom in 
the vector coefficients of the Fourier series minus AN + 1 
constraints. 

In the case of Nambu-Goto strings in an arbitrary num- 
ber of dimensions d one would span both right and left 
moving excitations according to 

1 

*' N (u)=l[E i R h2 (u)E?4 (44) 

i=N 

and 

1 

b' N (v) - J] E.R^Ejh'^ (45) 

i=N 

with the Ei given by the a choice of ( f43| ) appropriate to 
the desired number of dimensions. 



V. CONCLUSIONS 

We have generalised the solution to the unit magnitude 
constraint presented in |?J from three to four dimensions, 
casting it somewhat differently, in an effort to arrive at a 
general parametrisation of chiral superconducting strings 
with a finite number of harmonics. We have further 
shown how to construct loop solutions that satisfy the 
center of mass constraint and exclude overall orientation 
freedom. This result is useful because in studies of the 
properties of chiral loops, such as self-intersection and 
gravitational radiation properties, overall orientation of 
the loop is unimportant. 

Studies of chiral cosmic string loops with constant cur- 
rents [^1 and simple varying currents Q have been per- 
formed. Generally, however, we expect the current to 
be arbitrarily varying when loops are formed by intersec- 
tions involving different strings or if different segments 
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of the loop or string were at some point in causally dis- 
connected regions. This is a fairly generic situation and 
a study of the properties of more general chiral loops 
should account for these variations. 

Along the way, we have found that our modification 
of the method lends itself readily to a generalisation to 
arbitrary dimensions. We use such a generalisation to 
present solutions that could be useful in the investiga- 
tion of classical relativistic strings in higher dimensions 
as well as strings in 3 + 1 Minkowski space with currents 
and charges induced by Kaluza-Klein compactification 
p| when the back-reaction from the gauge fields can be 
considered negligible. 
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